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Abstract 

The intrinsic decay rate of orthopositronium formed in Si02 powder 
is measured using the direct 27 correction method such that the time 
dependence of the pick-off annihilation rate is precisely determined. The 
decay rate of orthopositronium is found to be 7.0396 ± 0.0012(stat.) ± 
0.0011(sj/s.)/is~^, which is consistent with our previous measurements 
with about twice the accuracy. Results agree weU with the O(a^) QED 
prediction, and also with a result reported very recently using nanoporous 
film. 
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1 History: Orthopositronium lifetime puzzle 



Positronium (Ps), the bound state of an electron and a positron, is a purely 
leptonic system, and the triplet (I'^-S'i) state of Ps, orthopositronium(o-Ps), 
decays slowly into three photons. Precise measurement of this decay rate gives 
us direct information about quantum electrodynamics(QED) in bound state. 

Three precision measurements [1, 2, 3] of the o-PS decay rate were performed 
at Ann Arbor, which reported decay rate values much larger, i.e., 5.2 9.1 ex- 
perimental standard deviations, than a QED prediction[4] (7.039934(10) ^s~^) 
corrected up to 0{a^). This discrepancy has been referred as 'orthopositronium 
lifetime puzzle', and was long-standing problem. To elucidate discrepancies, a 
variety of experiments have since been carried out to search for the exotic decay 
mode of o-Ps, resulting in no evidence so far [5, 6, 7, 8, 9, 10]. 

As some fraction of o-Ps inevitably results in 'pick-off' annihilations due 
to collisions with atomic electrons of the target material, the observed o-Ps 
decay rate Xobs is a sum of the intrinsic o-Ps decay rate Ao-Ps and the pick-off 
annihilation rate into 27's, Xpick, i.e., 

Kbs{t) = X:i., + Xp^ck{t)- (1) 

Xpick{t) is proportional to the rate of o-Ps collisions with the target ma- 
terials, i.e.; Xpick = n(Jav{t), where n is product of the density of the target, 
da the annihilation cross-section, and v{t) the time dependent velocity of o-Ps. 
Due to the thermalization process of o-Ps, this necessitates expressing Xpick 
as a function of time whose properties are dependent on the surrounding ma- 
terials. Thermalization process should be carefully treated civen in the cavity 
experiment [3]. Although pickoff correction is small in cavities, disappearance of 
o-Ps through the cavity entrance aperture has large contribution to Xobs- This 
disappearance rate is also proportional to v{t), as the same reason. Since the 
rate of elastic collision is extremely small in cavities, it takes much time, longer 
than 1 /xs, to thermalize well, and the disappearance rate still depends strongly 
on time. 

In previous measurements[l, 2, 3], Xobs^s were measured by varying the den- 
sities of the target materials, size of the cavities and also the entrance aperture of 
the cavities. The extrapolation to zero density or aperture was expected to yield 
the decay rate in a vacuum. As-,,, under the assumption of quick thermalization 
(shorter than 170-180 nsec) with constant o-Ps velocity. However, this assump- 
tion contains a serious systematic error as pointed out in reference[ll, 12]. 

We have proposed the following entirely new method [12], which is free from 
above-mentioned systematic error. The energy distribution of photons from 
the 3-body decay is continuous below the steep edge at 511 keV, whereas the 
pick-off annihilation is 2-body which produces a 511 keV monochromatic peak. 
Energy and timing information are simultaneously measured with high-energy 
resolution germanium detectors such that Xpick{t) / X^-y can be determined from 
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the energy spectrum of the emitted photon. Once a precise thermalization 
function is obtained, Xpick{t) will contain all information about the process. 
The population of o-Ps at time t, N{t) can be expressed as 

N{t) = AT^ exp (^-Xs,J^(^l + dt'^ . (2) 

Providing the ratio is determined as a function of time, the intrinsic decay rate 
of o-Ps, Aa-y, can be directly obtained by fitting the observed time spectrum. 

We obtained decay rate of 7.0398(29) and 7.0399(25) fis'^ independently [12, 
13], which are consistent with the non-relativistic QED calculation [4], and quite 
differ from the results obtained at Ann Arbor[l, 2, 3], 7.0482(16) 7.0516(13)^*3^1 
The observed Xpick (t) indicates that o-Ps thermalization is slow and it is serious 
systematic problem in all experiments using an extrapolation. In 1998, they rec- 
ognized that the incomplete thermalization makes the serious problem in their 
results[15], but they did not update results, and the discrepancy still remaining. 

There are also still several problems remained in our result[12]: (i) accuracy 
was 350 ppm, being worse than those of the other experiments [1, 2, 3], (ii) There 
was unknown systematic uncertainties before t start = 200 ns, and therefore to 
remove this uncertainty, final results were obtained using data after 220 ns, and 
(iii) systematic error regarding the Stark effect was not estimated. Improving 
results by considering these problems is described [14] here, and we report the 
final result using our new method. 

2 Experiment 

Figure 1 shows a diagram of the experimental setup; A ^^Na positron source 
(dia., 2 mm) with approximate strength of 0.03/iCi, being sandwiched between 
two sheets of plastic scintillators and held by a cone made of aluminized mylar. 
The cone was situated at the center of a cylindrical 50 mm-diameter vacuum 
container made of 1 mm-thick glass, being filled with Si02 powder and evacu- 
ated down to 5 X 10~^ Torr. Two different types of Si02 powder were used as 
listed in Table 1, with the biggest difference being in the mean distance between 
grains such that different pick-off ratio would be obtained. Using these powders, 
two 6-month runs were performed. 

Most of the emitted positrons pass through the scintillator, transmitting a 
light pulse to a photomultiplier(Trigger PMT) and forming Ps when stopped 
in the silica powder. Two high-purity coaxial germanium detectors (Ortec 
GEM38195) precisely measured the thermalization process. Energy resolutions 
were measured using several line 7 sources, with typical resultant values of 0.53 
and 0.64 keV at 514 keV. Lead sheets 2.0-mrn-thick were placed in front of 
each detector for suppressing contributions from simultaneous low-energy 7 hits 
from the 37 decay of o-Ps. Four large cylindrical Nal(Tl) scintillators (Scionix 
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Figure 1: Schematic diagram of experimental setup. 



RUNl RUN2 

primary grain size (nm) 7 7 

surface area {m^/g) 300 ± 30 260 ± 30 

density {g/cm^) 0.03 0.05 

mean distance between grains (nm) 340 200 

surface hydrophile hydrophobe 

Table 1: Characteristics of Si02 powders used in the measurements. 



127A127/5; 127 mm {(p) x 127 mm [t)) simultaneously measured the time and 
energy information from each decay. 

A new time-to-digital converter (TDC) was developed and employed for the 
present measurement. This direct clock-counting type TDC(2GHz) has a time- 
resolution of 0.5 nsec with known accuracy of 1-ppm. The time range for each 
channel is 32 fis, and the integral non-linearity, which is an important source of 
systematic errors, is controlled to be extremely small at < 15 ppm. To provide a 
systematic check, we used a 200-MHz internal clock-based TDC with 5 ns time 
resolution. 

The trigger PMT signal is fed into a fast leading discriminator whose output 
provides common start signals for the TDCs. Each output signal from the 
photo-multipliers to the NaI(Tl)s is fed into three ADCs and a fast loading edge 
discriminator that provides stop signals for the TDCs and gate timing for the 
ADCs. One ADC with a 3-fis gate width, measures the whole charge for the 
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duration of the signal. The other two ADCs eliminate pile-up events at the tail 
of the signal and base-line fluctuations, respectively. One ADC measures the 
signal charge with a narrower gate width (250 ns) comparable to the intrinsic 
decay time constant of the Nal(Tl) scintillator, while another ADC measures 
the base-line condition of the signal (180 ns width) just prior to the event. All 
gates for these ADCs are individual and gate timings are optimized for signal 
timings. 

Signal outputs from the Ge detectors are used for timing measurements and 
obtaining precise energy information. One signal is fed into a fast-filter ampli- 
fier (FFA) whose output is used as the stop signal for the TDCs, two auxiliary 
ADC signals, and three discriminaters whose thresholds are set at —50, —100, 
and —150 mV. Each FFA output is fed into another high-resolution TDC having 
fine resolution of 250 ps. Signal timings are determined utilizing signal shapes 
calculated from different threshold discriminating times and extrapolating to 
the intrinsic timing[14]. Good time resolution of 4 ns is obtained, as is efficient 
rejection power for slow rise signal components known to disturb Ge timing 
spectroscopy. Similar to the Nal(Tl) detectors, the two auxiliary ADCs effec- 
tively reject pile-up events by measuring the signal charge with a narrower gate 
and earlier timing gate. Precise time and energy resolutions can be obtained 
with this rejection scheme. The other output from each Ge detector is indepen- 
dently fed into a spectroscopy amplifier. Amplified with a 6-/is time constant, 
output is provided to a Wilkinson type peak- holding ADC that provides a 
precise energy spectrum. 



3 Analysis and discussion 

The ratio Xpick{t) / Xsj is determined using the energy spectrum measured by 
the Ge detectors. The energy spectrum of the o-Ps decay sample, referred to 
as the o-Ps spectrum, is obtained by subtracting accidental contributions from 
the measured spectrum. 

The 37-decay continuum spectrum is calculated using Monte Carlo simula- 
tion in which the geometry and various material distributions are reproduced 
in detail. For every simulated event, three photons are generated according to 
an order-a-corrected energy spectrum[20]. Successive photoelectric, Compton, 
or Rayleigh scattering interactions of every photon are then followed through 
the materials until all photon energy is either deposited or escapes from the 
detectors. The response function of the detectors is determined based on the 
measured spectrum of monochromatic 7-rays emitted from ^^^Eu, ^^Sr, and 
^^^Cs, with this function being used in the simulation. We refer to the obtained 
spectrum as the -spectrum which is normalized to the o-Ps spectrum with the 
ratio of event numbers within the region (480-505 keV). 

Figure 2(a) shows good agreement between the o-Ps spectrum and 87 spec- 
trum below 508 keV, where the pick-off annihilation peak is evident at the edge 



4 



of the 37-decay continuum. The lead sheets in front of the Ge detectors effec- 
tively suppress the contribution from events in which two low-energy 7's emit- 
ted from a 87- decay simultaneously hit one detector and that the 37-spectrum 
well reproduces such simultaneous events. Figure 2(b) shows an enlarged view 
of the observed o-Ps spectrum after subtracting the 37-spectrum, where good 
agreement with a detector response function is present. The ccntroid of this 
spectrum is 510.997 to.oH (keV), which is consistent with 511.0 (keV). This in- 
dicates successful subtraction of the 87 contribution, and the resultant peak can 
be regarded as pure pick-off annihilation samples. Obtained ratio of Xpick/^s-y 
is 0.01049(8) for time window of 150 - 700 ns. 

These calculations of the Xpick/^a-y are performed for various time windows, 
and its time dependence is shown in Fig. 3. Since the fractional energy loss of 
o-Ps per collision with Si02 powder and the collision rate are both dependent 
on its energy, the time dependence of the average kinetic energy of o-Ps at time 
t, E{t) can be derived from the Boltzmann equation 

^W) = -\/2mp,W) (W) - IkBT^ E ' (3) 

where nips is the mass of o-Ps, T room temperature, and ks the Boltzmann 
constant. The last term, the momentum transfer cross-section of Si02 is ex- 
panded in terms of E{t), i.e., the coefficients aj represent the effect of effective 
mass at the surface of the Si02 grain and mean distance between the grains. 
Since the pick-off rate is proportional to the average velocity of the o-Ps, the 
ratio 6{t) = Xpickit)/^^^ can be expressed by a differential equation[16], i.e., 

leit) = -c{9itf-9l)eitf^, (4) 

where C is a constant, 9ao = d{t 00), and the last summation term in 
Eq. 3 is replaced with an arbitrary real number /?. The measured Xpick{t) / Xs-y's 
are fit with this equation and figure 3 shows best fit results using the MINUIT 
library, where the pick-off rate cannot be assumed as constant, even in Si02 
powder where the collision rate is expected to be higher. 

Figure 4 shows the time spectrum of Nal(Tl) scintillators for RUN 1 with an 
energy window of 370 — 440 keV, in which a sharp peak of the prompt annihi- 
lation is followed by the exponential decay curve of o-Ps and then the constant 
accidental spectrum. The o-Ps curve is widely observed over about 1.2 /xs cor- 
responding to about eight times the o-Ps lifetime, because of the weak positron 
source used in this experiment and good suppression of accidental contributions 
by selecting the 7-energy. Observation in the wide range is important to under- 
stand time-dependence of A. To effectively eliminate pile- up events, a base-line 
cut condition was applied, and to further reject pile-up events, events with small 
differences between two ADC values (wide and narrow gates) were selected. 
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Figure 2: (a) Energy spectrum of o-Ps decay 7's obtained by Ge detectors. Dots 
represent data points in a time window of 150 — 700 ns, and the solid line shows 
the 37-decay spectrum calculated by Monte Carlo simulation. Shaded area 
indicate simultaneous hits estimated by the simulation, (b) Pick-off spectrum 
obtained after subtracting the 87 contribution from the o-Ps spectrum. The 
solid line represents the fit result. 



We fit resultant time spectrum using the least square method, i.e., 



Nobs{t) = exp{-Rstopt) 



£37 



N{t) + C 



(5) 




Figure 3: The ratio Xpick{t) / ^a-y are plotted as a function of time. Open circles 
are data points for RUNl and closed circles for RUN2. Solid lines represent 
best fit results obtained using Eq. 4. 



where Spick and 63-^ are the detection efficiencies for pick-off annihilations and 
37 decays, and Rstop is an experimental random counting rate representing the 
fact that time interval measurement always accept the first 7 as a stop signal. 
Apicfc/As-y is about 1% due to the low-density of the Si02 powder, i.e., the ratio 
of error propagation to decay rate is suppressed by a factor of 100. 

Figure 5 shows obtained fitting as a function of fitting start time for both 
runs, where values arc stable with respect to fitting start time except those 
before 80 ns. Since the fitting x^'s for both runs rapidly increase before 80 ns 
due to the tail effect of prompt peak, values a,tt = 100 ns are taken as the final 
results. The reduced x^'s a.t t = 100 ns are 1.006 and 1.051 for RUN 1 and 
RUN 2, respectively. The fitted values of decay rate are extremely flat after 90 
nsec, and are confirmed to be independent of fitting start time^ . The obtained 
decay rates are X^^ = 7.03991 ± 0.0017{stat.) ^5"^ for RUN 1 and 7.03935 ± 
0.0017{stat.) /is^^ for RUN 2, which are consistent with each other. As shown 

^The systematic increase in decay rate before t start = 200 ns observed in previous 
measurement[12] was eliminated using the new system. 
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Figure 4: Time spectrum of Nal(Tl) scintillators for RUN 1 within an energy 
window of 370 — 440 keV. The inset shows an enlarged view of the prompt peak, 
where good time resolution of cr = 2.2 ns is obtained. 

in Fig. 3, the thcrmalization process and pickofF ratio are different within each 
run, the good correspondence between two runs is obtained. This indicates that 
our method correctly takes into account thcrmalization and pickoff correction. 

Table 2 summarizes estimates of various systematic errors, and the details are 
shown in reference [14]. We just mention here about errors related to simulation 
and Stark effect. 

Monte Carlo simulation related errors: the predominant contribution to total 
systematic error is produced by uncertain normalization. That is, the number 
of pick-off events are determined by subtracting the normalized 3-7 spectrum of 
Monte Carlo simulation from the o-Ps spectrum, where changing the normaliza- 
tion factor affects the Xpick (t) / X^-y values and eventually propagates to the final 
result. Since the sharp fall-off of the 37-spectrum at 511 keV is solely produced 
by the good Ge energy resolution of a = 0.5 keV, this subtraction only affects 
the lower side of the pick-off spectrum such that improper subtraction results 
in asymmetry of the pick-off spectrum shape. Comparison of the asymmetries 
of the pick-off peak shape and the prompt peak annihilation spectrum is a good 
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Figure 5: Decay rates as a function of fitting start time. Small horizontal lines 
on the error bars represent the size of statistical errors solely due to fitting, 
while vertically extending bars include propagated errors from Xpick (t) / Xs-y de- 
termination. Open and closed data points indicate values for RUNl and RUN2, 
respectively. Dashed and solid lines show one standard deviation obtained at 
t = 100 ns for RUNl and RUN2, respectively. 



parameter for estimating this systematic error. The 1 a error is assessed as 
±99 ppm. Inhomogcncity of the powder affects on the detection efficiency. Si02 
powder density in the MC simulation is conservatively changed by ±10% al- 
though the uniformity is known to be within a few %; a change resulting in an 
error of < ±55 ppm. 

The Stark shift stretches the lifetime of Ps atoms, i.e., a perturbativc calcula- 
tion shows that the shift is proportional to a square of the effective electric field 
£; such that AAs^/Aa^ = 248-{E/Eof, where Ea = mle> /Ti" w 5.14x10^ V/cm. 
E is defined as the root-mean-square electric field sensed by o-Ps during its life- 
time. Calculations have estimated two contributions exist based on measure- 
ments of the electrical charge-up on the primary grains of silica powders and 
electrical dipole moment on the surface of grains[14]. The charge-up is partly 
intrinsic depending on powder specifications and partly due to positron deposi- 
tions from the /?+ source. The effect, however, is negligible in both cases, i.e.. 



Source of Contributions 


RUNl (ppm) 


RUN2 (ppm) 


TDC module dependence 






- Calibration 


< 1 


< 1 


- Stability 


2~3 


2 ~ 3 


- Integral Non Linearity 


< 15 


< 15 


- Differential Non Linearity 


Negligible 


Negligible 


Cut condition dependence 






- Base Line Selection 


-17 and +89 


-11 and +23 


- WD-NW condition 


-6 and +45 


-20 


Monte Carlo dependence 






Normalization 


±99 


±113 


- Relative efficiency 






of Nal(Tl) scintillator 


+7 


+7 


- Inhomogcncity of Si02 powder 


< ±55 


< ±88 


Other Sources 






- Zeeman effect 


-5 


-5 


- Three-photon annihilation 


-26 


-42 


- Stark effect 


-21 


-6 


- n = 2 excited state 


+19 


+19 


Total 


-120 and +153 


-151 and +148 



Table 2: Summary of systematic errors. See text for details of each item. 



on the level of 10~^ ppm at most. Silanol functional groups on the surface of the 
powder grain behave as an electrical dipole moment creating an effective field 
around the grains. Average densities are measured to be 2.5/nm^ and 0.44/nm^ 
for RUN 1 and 2, respectively. Accordingly, the effective field can be analytically 
calculated such that the contribution to the o-Ps decay rate is determined to be 
—21 ppm for RUNl and —6 ppm for RUN2. These estimations were confirmed 
using results from precise hyper-fine-structure (HFS) interval measurements of 
ground state Ps in silica powder[18], where the interval is proportional to the 
size of Stark effect. Considering the difference in powder densities used, the 
HFS results are consistent with our estimation. 

Other sources of systematic errors: Error contribution due to the Zeeman ef- 
fect is estimated using the measured absolute magnetic field around the positro- 
nium assembly (—5 ppm). Since the 8-7 pick-off process can only occur at a 
certain ratio, the calculated relative frequency (T37/CT27 ^ 1/378 is consistent 
with previous measurements [17], being —26 ppm for RUNl and —42 ppm for 
RUN2. The probability of the excited state (n=2) of Ps is about 3 x 10"'' [19], 
which could make the intrinsic decay rate 19 ppm smaller due to the low decay 
rate (gAs-y) 
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Figure 6: Historical plot of o-Ps decay rate measurements including present 
results. Vertical line shows the 0(a^ )-corrected NRQED prediction. Small 
vertical lines on the error bars indicate the size of errors solely due to statistics 
while error bars represent total ambiguities including systematic errors. 



The above discussed systematic errors are regarded as independent contri- 
butions such that the total systematic error can be calculated as their quadratic 
sum, resulting in —120 ppm, +153 ppm for RUNl and —151 ppm, +148 ppm 
for RUN2. 



4 Conclusions 

The decay rate of o-Ps formed in Si02 was measured using a direct 2j correction 
method in which the thermalization effect of o-Ps is accounted for and integrated 
into the time spectrum fitting procedure. Consistent results were obtained using 
two runs with different types of Si02 powders, i.e., and a weighted average gave 
Ao_Ps = 7.0396 ± 0.0012{stat.) ± 0.0011(s?/s.)^s"\ which agrees well with our 
previous result[12, 13], it disagrees with results obtained at Ann Arbor, i.e., by 
3.8-5.6 a [2, 3]. As illustrated in Fig. 6, our value agrees well with the NRQED 
prediction corrected up to O(a^) term[4]. We confirm our results obtained at 
1995 and 2000, and improve accuracy by factor of 1.8. 
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Just one week before this workshop, Ann Arbor reported new result[21] 
using nanoporous sihca film, which produces near-thermal energy o-Ps. They 
recognized the systematic problem due to unthcrmalizcd o-Ps (fast o-Ps) and 
suppress the kinetic energy of the emitted o-Ps using nanoporous silica film and 
positron beam energy. They obtain value of Aq-ps = 7.0404 ± 0.0010{stat.) ± 
0.0008(sj/.s. )/is^^, which completely agrees with our results as shown in Fig. 6. 
So we can conclude that the unthermalized o-Ps has contributed large systematic 
errors of about 1000 ppm and made discrepancy between the QED prediction 
as pointed out in references [11, 12]. Now two groups (two methods) obtain 
the consistent results and the orthopositronium lifetime puzzle is clearly solved 
here. The combined results of two different methods is 7.0400(10)/is~^, which is 
not enough accuracy(14Gppm) to examine the 0(a^ )-correction (210ppm). We 
need much effort to improve accuracy. 

We wish to thank Prof. S. Orito, T. Hyodo and Y. Nagashima (Univ. of 
Tokyo) for very useful suggestions for discussions. It is a pleasure to thank Prof. 
G.S. Adkins (Franklin & Marshall College) for calculating the 0{a) corrected 
energy spectrum. 
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